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Abstract 

A technique for evaluating the electromagnetic Casimir energy in situations involving spherical or circular 
£Nj ' boundaries is presented. Zeta function regularization is unambiguously used from the start and the properties 

of Bessel and related zeta functions are central. Nontrivial results concerning these functions are given. While 

■ part of their application agrees with previous knowledge, new results like the zeta-regularized electromagnetic 

■ Casimir energy for a circular wire are included. 

G\ 

1 Introduction 

Casimir effect problems have been showing a remarkably long-lived appeal since the day of their birth|l|, still 
stirring up intense activity all along the eighties to name just a few works), and reaching present-day 

topics (e.g. During all this time they have been object of many different approaches: stress-energy 



tensor |12[, Green function methods [|13|], multiple scattering expansions |l4j, heat-kernel series (| 15|-(l9[|), etc. 

In the present paper, we offer new calculations of the Casimir energy for an c.m. field in the presence of a 
sphere in D — 3, and of a circle in D = 2, with perfect-conductor conditions at their points. The contributions 
from inside and outside the boundary are separately studied. Further, in D = 3 we take apart the pieces associated 
to trasnverse electric (TE) and transverse magnetic (TM) modes in each case. Although this is lengthier than 
starting from the whole sum (since some pieces which cancel would then be eliminated at the beginning, and 
now we keep them until the end) we get as a bonus a useful decomposition into the four contributions. Thus, we 
find as a byproduct the effect coming from e.g. the interior part only or the Dirichlet modes by themselves. 

Unlike previous works on this subject involving zeta-functions at some stage or other (e.g. [^0[ 0), we 
adopt here what might be called a 'straight' or 'frontal' zeta- regularization approach right at the outset from the 
eigenmode summations themselves, as advocated e.g. in || 

We believe that our attitude is quite sound, as generic zeta function regularization pl| , |22] | is already a 
widespread technique. This method, often applied when boundary conditions affect a given quantum system 
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— like quantum billiards |2^, |24j — or field theory, is also adequate for dealing with finite temperature or curved 
spaces. Moreover, its variants have suceeded in problems involving nontrivial topologies or boundaries, like 
Kaluza-Klein style models [E5l Eq], or systems under the influence of external fields |27f or in interaction with 



material bodies (e.g. superconductors 28 ) — other applications are shown in |29| . 

Many of these advances came true by means of a precise knowledge of the spectral zeta function associated 
to each system, which is the key to the derivation of many results. Somehow, this happens also to the present 
paper: our calculation is indebted to the investigation started in |}0| about Bessel zeta functions^, |34|, and 
taken further in |j38f for the case of a purely scalar field inside a sphere, under Dirichlet boundary conditions. 
That work is here extended so as to include all the field modes appearing in the problems under discussion. It is 
also worth to note that a similar technique has been applied in Quantum Cosmology to study the Hartle-Hawking 
wave function of the universe [ f39| . 

After a brief survey on electromagnetism (sect. 2) general considerations about the physical problem of the 
Maxwell modes, and its spectral zeta function are made in sect. 3. First, we study in some detail the part 
associated to the internal transverse electric (TE) modes with particular emphasis on the zeta function for the 
zeros of the J v Bessel function and on the construction of the corresponding complete zeta function. Then, the 
same procedure is applied to the internal transverse magnetic (TM) part and its analogous zeta function. Next, 
the method is repeated for external modes in sect. 4. Comments about the result in these and related situations 
are made in sect. 5. Essential mathematical material concerning the derivation of partial wave zeta functions 
for external modes has been packaged into app. A. App. B contains an example of complete zeta function 
alternative calculation by using the coefficients of the heat-kernel series for the Laplacian. A comparison with 
previous estimates for D = 2 is made in app. C. 



2 Electromagnetism in D = 3, 2 

2.1 Neutral and perfectly conducting sphere in D = 3 

We shall briefly sketch the classical problem of an electromagnetic (e.m.) field kept within a cavity resonator 
bounded by a perfectly conducting spherical shell of radius a. The adequate conditions on the surface are 
n ■ B\ r=a = ft X E\ r=a = (where ft is the normal vector), in addition to the requirement of regularity in 
the interior. Then, the spherical e.m. waves which result from solving Maxwell's equations have radial parts 
(x r 1_D / 2 J„(fl ;)(wr) 1 where 

u(D,t) =1 + ^-1 (2.1) 

denotes the Bessel index for angular monentum I in D = 3 space dimensions. In the present cirsumstances, the 
solutions are divided into 'transverse electric' (TE) and 'transverse magnetic' (TM) ones (see e.g. Their 
possible associated frequencies (uj) are then determined by conditions on the radial parts which read, in each 
case, 



r l ~ D,2 J v ^i){ur) = 0, for region I TE-modes, (2.2) 



d_ 

dr 



D-2 l-D/2 t / \\ 

r r 1 J v { D j)(ujr) j 



= 0, for region I TM-modes, (2.3) 
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Clearly, when a particular solution is of TE type, toa has to be a nonvanishng zero of a Bessel function J v m 



for some /. (2.3) can also be written as 

(D/2 - l)J vm) (k) + U' HD>l) (k) = 0, k = wa, 



(2.4) 



which is a Robin (or standard homogeneous) condition with relative coefficients (D/2 — 1, 1). If one keeps in 
mind the Maxwell equations, all this is valid for D = 3 only but, when regarded as just two massless scalar fields 
obeying the Klein-Gordon equation and subject to Dirichlet and (a specific) Robin b.c, it holds for general 
The same boundary conditions apply to the external (region II) solutions and we shall come back later to this 
question. 

What is more, taking advantage of the D — 3 duality E <-> B between e.m. and colour gauge fields one realizes 
that, up to a global factor equal to the number of SU(A r c ) degrees of freedom, the present set-up is equivalent to 
a bag model for linear QCD0 including only internal gluon modes. 



2.2 Neutral and perfectly conducting circular wire in D = 2 

Two-dimensional spaces can be of interest since they are the scenario for systems such as quantum billiards p^, 
p3| , [24j] or anyons p5fl . In particular, the electromagnetic-Chern-Simons Casimir effect in a (2 + l)-dimensional 
spacetime has already been studied for parallel conducting lines in |3(| and for a circle in |37[] . What is explained 
in this subsection might also be obtained by the formalism set up in sects I and II of |37| when the Chern-Simmons 
part is absent. The components of the e.m. tensor in a (2 + l)-dimensional spacetime are specified in terms of 
its electric and magnetic fields 

/ Ei E 2 \ 



± flis 



-E x 
-E 2 



B 
-B J 

In D = 2 there is room for only one B-component. Actually the magnetic field has become a scalar rather than 
a vector. We impose the Maxwell equations in the vacuum (absence of charges and current density) 



dpF*" = 



d a F* a = 0, 



where F* a = ^^"F^ is the dual tensor to the e.m. field. More explicitly we have F* = B F* 1 = 
F* 2 = Ex. From <U) we get 

rot E = B. 



where we understand that rot E = d\E 2 — d 2 E\. From equation (2.5) we draw 



div E = 
E x + d 2 B = 
E 2 - diB = 0. 



(2.5) 

(2.6) 
-E 2 and 

(2.7) 

(2.8) 
(2.9) 
(2.10) 



1 The possible relevance of the space dimension as a perturbation parameter was noted in Jjij . 

2 by 'linear' we mean that the nonlinear f ahc A^AY term is omitted from the field-strength tensor, as usual in such approache 
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(2.7)-(2.10) are the complete set of fundamental equations. Taking the time-derivative of 82 of (2.9), 8\ of 

(2.10) and combining the results, one finds V 2 B — B = .i.e. 



□5 = 



(2.11) 



which shows that B obeys a Klein-Gordon equation. Now, we give the form of these expressions when the fields 



have a time dependence of the form E(x, t) = e(x)e lut and B(x, t) = b(x)e 



rot e 
div e 
iuei — 82b 
%LOt2 + d\b 



-iu/b, 



0. 

0, 

0. 



(2.12) 
(2.13) 
(2.14) 
(2.15) 



From this set of equations — or straighforwardly from (2.11) and the time dependence of B(x,t) — we get that 
(A + uj 2 )b = 0. Once this equation is solved, the electric field is given by equations (2.14), (2.15). The ensuing 



solution may be easily seen to automatically satisfy expression ( 2.12| ) . Now we arrive at the question of boundary 
conditions. Ancient lore tells us that the electric field must be orthogonal to the surface of a perfect conductor; 



if the normal vector is given by (711,712), this condition takes the form = ri\E2 — ^2-^1, which, using (2.14) 



(2.15) is equivalent to = -&B. To sum up: the problem is reduced to that of a scalar field B which satisfies 



the typical Helmholtz equation and satisfies Neumann boundary conditions. In our notations, for D = 2 one has 
a(2) = echoing the conversion of the Robin b.c. into a purely Neumann condition. 

The reader may check that all the modes generate a non-zero total charge on the boundary (which is given 
by the scalar product of the electric field and the normal vector to the boundary) . 



3 Internal modes 

Within these types of set-up, zero-point energies emerge from vacuum mode-sums of the type Tic — ^""^ ui n , and 

n 

give rise to the celebrated Casimir effect jl], ||. Note that the summation extends over all the w„'s in the set of 
eigenmodes. As a result, sums of this sort do usually diverge and call for some regularization to make sense of 
them. 

At this point, we introduce the usual spectral zeta functions, which we denote by 

o<( S ) = 5>r, c^w = e(t) S ' (3 - 1} 

\x is an arbitrary scale with mass dimensions, often used to work with dimcnsionless objects. As they stand, these 
identities hold only for Re s > so, being so a positive value given by the rightmost pole of (m(s)- However, such 
a function admits analytic continuation to other values of s, in particular, to negative reals. 

The finite part of the vacuum energy — Eq — can be found by combining zeta-regularization of the mode-sum 
and a principal part prescription from ref. |(J . Following that work, one may put 

"1 



E c (fi) = PP S 



7ic/uCm(s) 



(3.2) 



where PP denotes principal part. One should be aware that the whole, observable, physical energy includes other 
terms which have to do with the couplings of the bag Lagrangian — see the discussion on this point in the same 
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reference. (From here on, we adopt the typical QFT units which make H = c = 1). Obviously, for this procedure 
to work we must know how to obtain the analytic continuation of Cm( s ) to — at least — a small part of the 
negative real axis. 

In order to proceed, we shall introduce 'partial- wave' zcta functions for fixed values of the Bessel index v. We 
define the zeta function for the zeros of J„ as (see also [Sol FHl ) . 



Cl' V (s) = J2jZi far Be a >1, 



(3.3) 



where j vn denotes the nth nonvanishing zero of J„. The 1,2? label has been added as a reminder that this 
comes from eigenmodes in region I dictated by Dirichlet-type b.c^] (Discrete versions of the Bessel problem, their 
solutions and associated zeta functions have also received some attention in p3[). Analogously, let 



6 n (s) = J2 k ^ , for Re s > 1, 



(3.4) 



with k vn denoting the nth solution of eq. ( |2.4[ ) for a given v. 

Reconsidering the same problem in D-dimensional space, taking into account the degeneracy d(D,l) of each 
angular mode in D dimensions, we define the 'complete' spherical zeta functions 



Cm 0) 



E w)E^(W» = flS E *(ao£5, («), 



n=l 



(3.5) 



l=l mia n=l l=lmin 

Imin is the minimum value of / and, for gauge fields in D — 3, l m i n = 1. If we consider scalar fields, then / m ; n = 0. 
As for d(D, I), we find its value in [Q and put 

_ fcmax(-D) 



d(DJ) = (21 + D -2) 



(l + D-3)\ _ \ 
l\(D- 2)\ ~ (D-2)! 



Z3-2-2fe 



fc=0 



where 



D — 3 

— - — for odd L> > 3, 

fcmax(2?) = ^ 

— - 2 for odd D > 4, 
and the form of the Ak (D) 's can be read off from E{| (see also pq] ) . 



(3.6) 



(3.7) 



3.1 Internal Dirichlet (TE) modes 
3.1.1 'Partial-wave' zeta function 



By (3.2), computing the Casimir energy calls for the knowledge of the Bessel zeta functions (3.3) at s = —1, while 



the complex domain where (3.3) holds is bounded by Re s = 1. Fortunately for us, Cl'^i 3 ) admits an analytic 
continuation to other values of s. What is more, in refs. |3(J and |3^| we showed how to obtain an integral 
representation valid for — 1 < Re s < 0, which reads 



.7 7"), , S . ITS 
7T Z 



dx x 



In 



/ 27rxe x Iv(x) 



, for -1 < Re s < 0. 



(3.8) 



In the mathematical literature, this object taken at even integer s is sometimes referred to as the Rayleigh function |33| 
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Yet, we shall have to work out ( p.8[ ) in order to obtain an equivalent representation more amenable to numerical 
calculation. The first step is to rescale x —* vx. Afterwards, we will perform a subtraction procedure on the 
resulting expression. The aim of that is the reduction of (s) to some elementary functions of s plus an integral, 
whose integrand should be relatively easy to express in terms of uniform asymptotic expansions^ (u.a.e.'s). 



The piece which we will subtract and add to the integrand of (3.£) is 



In 



(l + x 2 ) 1 /* 



1,27 i (1 + X 2 ) 1 / 4 t r, 

(J 1 m -= h a 2 v{rj{x) — x) 



<j{ = -1, 



0+ = +1- 



where f](x) is the function appearing in the known u.a.e. of the Bessel function (see e.g.|46|): 

x 



[x) = \/l + x 2 + In 



V 



and a\'' D , <j^° have been introduced for future convenience. 



1 + Vl + x 2 ' 



The added term is then separately integrated using the intermediate steps 
dx x~ 



dxx s (v( x ) ~ x ) 



4s sin ^ ' 



dx x 



In- 



1 

2? 



s + 1 
2 ' 2 



1+ Vi+^ 



+ / dx x s l (s/l+x 2 - x) 



B 



s + 3 



-,s 



which lead to 



„ S . ITS 

-v — sm — 

7T 2 



1,27 



B 



s + 1 s 
2 ,_ 2 
s + 3 



s + 1 



/>OC 

j o 



r(x)r(w) 

where, as usual B(x,y) = — ; -, and 

T(x + y) 



£ I,T, (v, x) ee y/2nv{\ +x 2 f l4 e- v ^ x) I v (vx). 
The advantage of this new representation is that the u.a.e. of the In argument has reduced to 



k=l 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



where the it^'s are known polynomials listed in books like |4q ]. Expression (3.12) is also handy by the way in 
which the pole at s = — 1 is exhibited. The singularity at this point is caused by: 
i) the B functions with one argument equal to 



ii) the large- a; behaviour of the integrand; by (3.14) one has 



In x)) = -t(x) + 0{t 2 {x)) = — + 0[- 



8vx 



*also called Debye expansions 



<» 



which asymptotically yields a logarithmic divergence. In fact, after integrating it has the same look as i), since 



dx : 



t{x) m = -B 



(3.15) 



in this case with m = 1. 



The calculation of the integral in ( 3.12 ) can be mentally divided into two stages. First: necessarily one has to 
delete and separately add the part responsible for the divergence in ii). Second: to make things numerically easier, 
and eventually extract some infinities which will appear later in the 'complete' zeta function, it is convenient to 
keep subtracting and adding more terms of the expansion of ln[£ I,I '(i/, x)] 



Thus, taking several terms in the series, we write the In function in (3.12) as 



In [C l ' v (y,x)] ~ In 



E 

fc>i 



u k (t(x)) 



E 

n>l 



K D {t{x)) 



where the U\i v 's are given by 



U 
U 
U 

U A 



t 

I 

16 „ 
25 i 3 

1384" 
13 1 4 

128 



5t 6 

~24~' 
3t 4 



531 1 

~64C 
71 t h 

~32 



5t e 
16" 



221V 

f ^28" 

531 t g 

64 



1105 1 

~ 1152 
339 t™ 

32 



!) 



565 1 



12 



128 



(3.16) 



(3.17) 



From them, we form the quantities 



(3.18) 



The expressions for these J'^' D (s)''s are easily obtained from the U^ T '(t) , s in ( 3.17 ) by application of ( |3.15 ), as a 
result of which it is enough to make the replacement 



With this, we may write 



2 V 2 ' 2 



/ dxx— 1 ln[£^(u,x)] = S l f(s, v) + £ ^lJ£Z, 

•/ „ i 



(3.19) 



S^ p (s,i/) 



U^(t(x)) 



(3.20) 



the key point being that S l ! ^ D (s,i') is a finite integral at s = — 1. Of all the ^^(sj's, J L ' u {s) is special as it 
contains the only contribution to the s = — 1 pole coming from the integral in ( [3.12| ) (i.e. the outcome of the 
'first stage'; the rest is produced by the 'second stage'). For this reason we change the notation to 

8 1 W 24 I 2 ' 2 



r I,X>/_\ P*'" D / s + 1 s \ . ^I,X>/ \ LB 1 



-B 



(3.21) 
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Here the symbol p l ' V has been introduced to make easier the expressions of the forthcoming cases. Then, we 
have 



1 



-!>£>.. -s 



_.S . ITS 

-v — sin — 

7T 2 



I 1 o / S + 1 S 



5 



5 + 1 
2 ' 



+2 S - 1 S 



s + 3 



(3.22) 



l^B ('-±1,-1)1 
2 F I 2 ' 2 / 



1 - Av 2 1 



n=2 

0((s + 1)°), which gives the right value of the 



Laurent-expanding near s = — 1 we find Cv{s) 

07T S + 1 

residue at this point |34|, [33j|. Similarly, by Taylor expanding close to s = we get ( u (s) — — i f ^ + ^ J +0(s), 
i.e., although we started from a representation valid for — 1 < Re s < 0, the correct value of £„ at s = is also 
recovered. 



3.1.2 'Complete' zeta function 

Next, we go on to the D-dimensional problem. This is done by inserting both ([0]) and ( 3. 22] ) into the first 



formula of (3.5). Afterwards, we may trivially interchange the /-summation and the fc-summation, since the 
second is finite. Once we have done so, we are left with sums of the type 



£ i/(A0~* = C»U T -i + J 



D 



(3.23) 



(recall ( |2.l| )) denoting the Hurwitz zeta function. Thus, we arrive at 



Cm ( s > 



2a 



(D 



1 



-a^Of I -D + 2 + 2k + s, - - 1 + /„ 



- £ (-l)^P) 
fc=o 



S . 7TS I it, 



1 R ( S+1 



x( H -D + l + 2k + s, 



D 



1 + L 



J] sjf ( a> i/(u,0)i/(A0 



D-2-2fc-s 



p» s (i±i,_i) c „(_ D+3+2t+ 4 



1 ~h ^min 



fJi'^^Cn (--D + 3 + 2fc + s,^ 



1 + L 



E ^( S X" + 2 + 2fc + s + n,^-l + uj | 



(3.24) 



Examining again the origins of singularity at s = — 1, we find 

a) The ones already present for £„, visible now as ~ B (^rp, ■ ■ •)• 

b) New pole contributions when the first argument of any of the present Hurwitz zeta functions equals one, 
including the terms with n and k such that n = D — 2k in the last sum. 
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On the other hand, at s = —1, the series S]f(s, v(D, I)) v(D, l) D ~ 



2-1k-s 



appears to have a rather slow 



numerical convergence but, since its net contribution is actually little (and the larger N, the smaller it gets) in 
practice it suffices to compute until an accuracy of few digits is achieved. 

3.1.3 Three-dimensional space 



As a result of ( p.7|) , for D = 3 one has fc m ax(3) = and the fc-series in ( 3.24 ) reduces to the k — 0-term, with 
Aq = 1. Further, we consider the description of the electromagnetic TE modes and therefore set l m i n — 1 (for an 
ordinary scalar field Z m ; n = 0). We split Cm ( s ) i n t° two pieces: CJvu( s )j containing the part directly evaluable 
at s ~ —1, and Cm2( s ) which includes the 'mixed' terms with both singular and finite contributions. Then 



Ca?(«) = Cmi(s) + Cm2( s ) 



*I.X> 



(3.25) 



where 



and 



&)=2a s 



2a s sin ■ 



1 ( i 3 



--sm T U 2 2 B( — 



-s ) ( H [ s - 2, - 



+ £3 



l>i 



N 



1-s 



+ JI' V ( S )CH U- + ]T Jn' V (sK* (n + 8-1, 



71^2 

n^3 



s 
2 



B 



Ch [s - 2, • 



(3.26) 



(3.27) 



Next, we will compute Cjvu( s ) directly at s = —1 and Laurent-expand Cm2^ s ) near s = ~ 1- With the specific 



values of <t 1 1 ' V , a 1 ^ , p I:V from (3J5), ( 3.21 ) and taking N — 4 subtractions, we get 



C^(-i) 



719 



1053 7T ■- i 

5760 ' 8192 1 65536 



35 7T 3 „ t v ( , 1 



(3.28) 



where, after a numerical calculation (with AT = 4), we have found 



1=1 



-1,1 + 



0.00024 



(3.29) 



The part containing the singularities gives rise to the following Laurent expansion 



Cm 9 («) = - 



2 

315 \s 
18457 



1 



1 



■In. 



60480 20160 
+0(s + 1)] . 



229 11 , „ ( „ 3\ 1 , 

7 ha2 - Cjt -3, - + -C w 

672 Sff I ' 2/ 4 SH 



(3.30) 



Here the prime denotes derivative with respect to the first argument. To complete the desired numerical evalu- 
ation, we still need the values of (' H ^— 3, and (' H ^— 1, , which are found from the relation C,h (s, = 
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— 2 S + (2 S — l)Cfl(s) — where means the Riemann zeta function — and from the knowledge of Cr(~3) ~ 0.005378 
and Ck(-l) ^ -0.165421. 



Equipped with all this, we are able to obtain 



Cm ( s ) 



315tt 



l 



In a 



The residue of the pole at s = — 1 is Res 



t M D (s),s = -l 



0.27069 + 0(s 
2 



3157TO 



(3.31) 

, Bearing in mind that a 1/2-factor appears 



when going from the zeta function for the Maxwell eigenmodes to that for the Laplacian spectrum, we realize that 
this agrees with the heat-kernel expansion of that operator in the Dirichlet case (see e.g. PI). The existence 
of this pole indicates that the Casimir energy under the present conditions is still infinite after zeta funtion 
regularization, and its divergence cannot be removed until the application of the PP prescription, which in some 
sense amounts to renormalizing. ( The issue of infinities for the bag model was considered, from the cutoff 
viewpoint, in |47fl . A full discussion about the essence of divergences in Casimir energy problems was supplied in 

S ) 



Repeating the calculation for i m j n = one obtains 



>scal,I,Z> 



a 



1 



+ lna + 0.00889 + 0(s + l) 



(3.32) 



315?r \s + l 

which corresponds to the internal modes of a true scalar field. In fact, this value may also be found by adding to 



(3.31) the 'partial wave' contribution of the I = mode alone, which is 



1 



1 



~Cl%(-l) = ^Cfl(-l) = ~ * -0.26180^ 
a ' a 12a a 



(3.33) 



(Note here that, since J\/2( x ) °t x^ 1 / 2 sin a;, Ci'/2( s ) = 7I " _s Cr( s ))' As one can check, the figures match. Moreover, 
the importance of the lower-lying region of the spectrum is manifest, since the / = part is almost as large as 
the rest, but with opposite sign. An approximate calculation of ( 3.32| ) based on the heat-kernel expansion of the 
Laplacian is given in app. B. 



3.1.4 Two-dimensional space 

One can now think of the circular wire problem in the plane. Given that we are still dealing with a Dirichlet 
field, this part will not be needed for the D = 2 e.m. Casimir effect, and is done just for completeness. Now 
d(2,l) — 2 for I > and d(2, 0) = 1. With the methods in the preceding subsubsections, i.e. carefully taking 
(|3.24|) for D = 2 and Z min = 1, we find: 



/-seal, I, T> / \ 1 

Cm "(-) = - 



16 



1 



In a + 0.02436 + 0(s + 1) 



(3.34) 



128?r Vs + 1 

Observe that i/(2, 1 = 0) = stops us from making the rescaling x — * vx and applying u.a.e.'s. A way of dealing 
with / = parts in D — 2 will be shown later. 



3.2 Internal Robin (TM) modes 



Here we outline the changes when considering the eigenmodes obeying the Robin b.c. ( |2.4| ). The analogue of 



(3.8) for these conditions is (see |3l| ) 



C^( S ) 



s ns 
— sin — 

7T 2 



dx 



In 



(xl' v (x) + al v (x)) 



for -1 < Re s < 



(3.35) 
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We shall take a(D) — D/2 — 1 in accordance with (2.4) and thus a(3) = 1/2. The calculation will be performed 
by a subtraction procedure similar to that for the TE modes but, now, the piece which we remove and add to 
the integrand is (instead of (|3.9|) ) 



-s-l 



hi 



(1 + Z 2 ) 1 / 4 ^(^j,)^) 



(l + a; a)i/4 LW 
(Tj In = h cr 2 ^(77(2;) — xj 



1,11 . -, 
fTi = +1, 



1,71 



= +1. 



(3.36) 



The ensuing expression is similar to (3.12) but for these changes 



a) 1 tie signs a 1 , er 2 are replaced with a{ , cr 2 . 



b) Instead of ( 3.13| ), for the present case one has 



C l ' n {u,x) = V2^- \r UI e-''^xI' v {yx) + a— -y^l + i 3 ) 1 ^')/,^) 

(1 + iz^) 1 ' 4 i/Vl + s 

Then In [£ I,K (i^, a;)] is expanded by taking advantage of the u.a.e.'s of both V v (yx) and I v (yx) (see again) 



where 



mfr^^a;)] - In 



U\' n {t) 



fc>i 



it>i 



E 



^(t(x)) 



3 \ 7t 3 
"sH** 24' 



3 a a 2 
16 + 2 ~ ~2~ 



t 2 + 



21 


3 a a 2 


a 3 ^ 


~128 + 


~8 2 


+ T, 


1463 i 9 






1152 ' 
27 


3 a a 2 


a 3 


~128 + 


~8 2 


+ T 


/ 733 
l _ "64 


Ala 
+ 8 " 


aA t 



5 a 
8 ~ 2 

3 i 



a 

T 



7t 6 



16 ' 

869 5 a a 2 
640 



— + — t° + 
4 2' 



109 23 a 3 a 2 
~32~ 8~~ ^ 



/441 21a 
l~32~ 8~~ 



2 

707 1 12 
128 ' 



315 7a 
128 + ~8~ 



(3.37) 



These formulas parallel eqs. (3.16) and (3.17). Looking at the linear t-term in the first of (3.37) we realize that 
the analogue of p l,v in the previous case (eq. ( [3.21 )) is 



IK ° 
P = a ~ o- 



(3.38) 



The expression for Cl ,7i {s) is like (3.22), replacing all the quantities with I, ^-superscripts by the corresponding 
ones obtained with I, "^-superscripts. When doing so, we must bear in mind that in the case under consideration 



dxx- 3 - 1 I In [£h n {v, x)] - 

{ n=l 

(which is by construction a finite integral) and 



Uk K (t(x)) 



^■i,Kt s 5 is + 3 s 



dxx-^U^itix)), n>2 



(3.39) 



(3.40) 



For studying the complete spherical problem, we must deal with the Ca^( s ) m the second line of (3.5), and apply 



the same formula ( |3.6] ) for the d(D,iys, obtaining a similar series of Hurwitz functions, integrals, etc. Later, we 
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specialize it to D — 3, (therefore a = 1/2) and Z m i n = 1- Before going on, we separate C, M ( s ) = Ctvii( s ) + Ca<2( s ) 
following the same finiteness criterion as in the decomposition (3.25). Then, the resulting Cmi( s )i Cm^( s ) 



are obtained from the Ca4i( s )j Cm2( s ) m formulas ( 3.26| ), ( |3.27 ) by just replacing all the objects having 1,2?- 
superscript by their counterparts with 1, 2C-superscript. We also need the numerical calculation 

2 



£3 

i=i 



1 



1 + 



1 



0.00012 



(3.41) 



(analogous to ( 3.29| ), and done for N — 4 too). After adding C^u( — 1) an d the Laurent-expansion of Cm^( s ) near 
s = — 1, wc find 



A^r \ 1 
Cm (s) = - 



-J- ( -J— +lna) - 0.10285 + 0(s + l) 
45?r V s + 1 / v ' 



(3.42) 



For a scalar field, one can either repeat the calculation with Z m j n = or separately add the contribution from 
this mode. The second option is easy because, using J\i<2(x) oc x -1 / 2 since, our Robin condition for v(3, 0) = 1/2 



just reads x 1 / 2 cosx 
1 



0. Its nonvanishing solutions are 7T (n + \ t 



0,1,2,. 



and therefore C^i/2( s ) 



tt ( s, — ] . At s = —1 one has 



a \ 2 



TT 1 

— = 0.13090-. 

24a a 



So, after adding this part, 



>scal,I,7?./ \ 1 

C M ■ ' I s ) = - 



3157T 



+ lna + 0.02805 + (9(s + l) 



(3.43) 



(3.44) 



The situation D = 2, a = may be similarly studied, arriving at the complete zeta function for the Neumann 



modes with L 



1: 



Cm ( s ) = C M ( s > 



a=0 



48- 5tt 
128tt 



In a 



0.17883 + 0(s + l) 



(3.45) 



4 External modes 



As already commented, the spectrum of modes — w's — is determined by the effect of the problem's conditions 
on the radial part of the wave solutions (in QFT language we would say 'the field solutions' ). This applies both 
to the solutions in the interior (region I) and to those outside the spherical surface (region II), which shall be 
now considered. There are at least two possible approaches, both of them leading to the same result: 
1. First, we sketch the simplest one. If we demand that the external solutions behave like outgoing radial waves 
~ e %ujr for r — > oo, (though the discussion might be repeated with asymptotically ingoing partial waves as well) 
their radial parts will be just oc r x ~ D l 2 H^ D ^(wr), with Hu denoting the first Hankel function (ingoing waves 
would just be the complex conjugate, thus involving H^> * = H^). Following this cue, we repeat the contour 
integration procedure of ref. f30fl , which gave the representations (^8), ( 3.35 ), but putting now H„ instead of 
J v at the outset. Since the calculations have the same nature, we do not go through! them here; perhaps the 
only point worthy of remark is that wherever in ref. ]3l| we took advantage of properties like 



argz < 7T, 



J u iz) <~ \ — cos z — 

" 7TZ ' 













(- 





















J„(e™ /2 z) = e W7r/2 /,(z), -tt < argz < |, 
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now we have to make use of 



0-("+5)i) 5 |z| > 1, -7T < argz < 2tt, 
2 

— ( 

7T 



HP(e*"' 2 z) = -de-^' 2 K v {z), -vr < argz < - 



(and their equivalents for conjugates and derivatives). Doing so, one obtains formulas (4.2) and (4.11) below. 
2. Another reasoning, physically more transparent, is to imagine a larger sphere of radius R, enclosing the one 
of radius a, on whose surface we impose conditions as well (in the siprit of ref.pl]]). Once their partial- wave 
zeta function has been constructed, the R — ► oo limit is obtained. Not only is the result i?-independent, but it 
coincides with the outcome of the method 1 as well. The full process is explained in app. A. 

Once the desired 'partial- wave' zeta functions Cu' V ( s )t Ci I,K ( s ) have been obtained by either method, taking 
the Bessel index v(D, I) and the degeneracy d(D, I) of each I in D dimensions, we construct the 'complete' 
spherical zeta functions 

oo 

= « s E d(Al)£(£ o 00, 

(4.i) 

i = imi„ 

When considering the e.m. field, we will eventually set Z m ;„ = 1 (for a scalar one, / m ; n = 0). 



4.1 Dirichlet (TE) external modes 

The adequate representation of the 'partial-wave' zeta function analogous to ( |3.S| ) for region II is 
e<%) = -sin^ I™ dxx-*- 1 In \l^e*K v { 

IT Z J Q V 7T 



>(*) 



for -1 < Re s < 0. 



(4.2) 



After rescaling x — > vx, a subtraction procedure similar to that for the internal modes will take place, but now 
the piece which we remove and add to the integrand is 



In 



As compared to (|3.9|), this amounts to a sign flip in v. Therefore 



(4.3) 



11,75 



11,73 -, 
(To = — 1. 



The ensuing expression is similar to ( [3.12 ) but with the above <7i, Oi instead of those, and the function £, I,T> (u, x) 
appearing there replaced now with 



C n ' v {u,x) = J— (l + x^e'^K^vx). 

V 7T 

In [£ ii,i '(j/, x)] is expanded by taking adavantage of the u.a.e. of K v {yx) (see [|6| again) 



\n[C n ' v {v,x)] ~ln 



i + E(-d 



k u k (t(x)) 



k>l 



E 

n>l 



ZC^(t(x)) 



where, by obvious parity reasoning with respect to the I, T> case (compare the above expression with (|I|)) 

U l ^ v {t) = (-l)X^(O- (4-4) 
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As a result of (O) 



Hence, by virtue of ( 3 . 18)) and of (4.4), 



J* l > v (s) = (-in»(s), n>2. 



(4.5) 



(4.6) 



With these elements, we can already construct the analogue of ( [3.22| ), say Ci 1 ' 23 ^)) f° r the external Dirichlet 
modes by 'superscript substitution'. The resulting expression contains 



S N ' (s,v) 



N 



dxx-*- 1 \ \n[L n ' v (v,x)} 



U^(t{x)) 



(4.7) 



which is a finite integral. When studying the complete spherical problem, we must handle the Cm^ X s ) m the 
first line of (4.1), and apply again formula (3.6) for the d(D,l)'s. Then, the emerging complete zeta function 
is like ( 3.24 ), but with all I, P-superscripts turned into II, D-superscripts. After setting D — 3, l m i n — 1, we 
Laurent-expand around s = —1, finding that everything is of the same sort as in the internal case, except for the 

new quantity j>2 S^' V (-1,1+ ^) (l + 7j) - -0.00054 (here calculated for N = 4 too). So, after gathering 



everything together one arrives at 

/•n,z>/ \ 
Cm ( s ) = 



315tt V s + 1 



1 



In a 



0.00326 + 0(s + 1) 



(4.8) 



Redoing the calculation for Z m i, 
/ = mode is not contributing in this case 
When D = 2, l m [ n = 1, we get 



0, corresponding to a scalar field, we obtain Cm ' l ^(s) = Ca^ I '( s )i i- e - the 



(s) = - 



16-7T 
1287T 



1 



+ lna + 0.00501 + 0(s + 1) 



(4.9) 



Now, the net contribution from the I = Dirichlet mode in D = 2 will be found by joining the inner and outer 



partial wave zeta functions (p.8|), (4.2). This sum has the effect of cancelling the s = — 1 divergences as a result 
of which we can numerically integrate, obtaining: 

1 



1 r 

a 



1 

air 



dx ln[— 2xIq(x)Kq(x)] 



-0.02802- 



Actually, by a calculation based on a slightly different representation of (j; v (s) valid for v = 
that the part of ( 4.1C| ) coming from the internal modes is 

1 1 



(4.10) 
we know 



1 J,x> / \ 1 
-C=o* = - 
a a 



-0.01451 + O(s + l) 



_8tts + 1 

4.2 Robin (TM) external modes 

Applying any of the above referred procedures we find the 'partial- wave' zeta function representation 



Cu' (s) = - sin — 

7T Z 



The subtracted part will be 



dxx s 1 In 



— e x (~xK'Jx) -aKJx)) 

TTX 



for -1 < Re s < 0. 



In 



(1 + x 2 ) 1 / 4 



-v (rj(x) —x) 



(4.11) 



(4.12) 
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ii, n 



-1. 



Thus, paralleling ( [3. 36 ), 

Oi ■ = +1, o 2 

With respect to the situation described as I, TZ, this amounts to a sign change in v and the corresponding parity 
reasonings apply everywhere ( as when going from 1,2? to 11,2?). In the present case we have to deal with 

£»•>,*) = -JE 1 e^xK' v {ux) - a^2==J^(l + i»)'^K 
V 7T (l + x z 1/4 z/\/l + a; 2 V 7r 



N 



dx : 



-s-l 



U^{t{x)) 



(4.13) 



n = l 



where Z#' w (t) = (-l) n Z# w (t) and therefore p 11 ^ = -p 1 ^, J^ n {s) = J^''*(s) = (-l) n ^' w («), 

n > 2. 

Then, we construct the complete zeta function, which we calculate for D = 3, a = 1/2. Taking into account 
that (for N = 4) ^ S* 1 ' K ^-1, Z + ( Z + ^) ~ -°- 00041 > 

we are able to write the Laurent expansion near 

s = -1: 

2/1 \ 

" ' La a -0.07223 + 0(s+l) 



. . . (4.14) 

45tt V s + 1 ' 

As for a possible Z = contribution when considering a scalar field we see, like in the external Dirichlet case, 
that this mode changes nothing, i.e. (^^'^'^(s) — Cjvf K ( s )- 
When D = 2, a = 0, 



1 



a=0 a 



48 + 5tt 

1287T 



1 



1 



In a - 0.03804 + 0(s + 1) 



(4.15) 



As in the D = 2 Dirichlet case, the total contribution from the I = mode alone follows from adding the inner 



and outer partial wave zeta functions which are now (3.35) and (4.11). Again, the s = —1 divergences cancel 
and we can find, by numerical integration: 



1 r 

a 



cS(-i) + S(-i) 



= — / dx ]n[-2xI' (x)K' (x) 
a =o an J 



-0.50704^ 

a 



(4.16) 



5 Discussion 



5.1 D = 3 

One can consider the contribution to the e.m. Casimir energy coming from the interior of the sphere only, i.e. 



the zero-point energy of a photon bag banning the existence of outer modes. Taking into account (3.2) and the 

■ ln(afi) + 0.08392 . 



sum of (3.31) plus (3.42) one finds 



E, 



c . m . I 



3157T 



(5.1) 



The logarithmic term, depending on p, has to be viewed as a remainder of the renormalization process implicit 
in the prescription adopted ||. The second article of refs. [^(| gives, for the Casimir energy due to the vector 
gauge bosons in the interior of such a bag 



■In- 



0.08984 



(5.2) 



315tt 8 

found by the 'energy method' of ]13| . 8 is a cutoff arising from the non-coincidence in time of field points, 
linkable to the nonzero 'skin depth' of a real — and not purely mathematical — surface. Incidentally, Riemann 
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zeta functions were also employed in part of the calculations in 1 20 , although the initial regularization approach 
was essentially different from ours. Therefore, altough they are within a 6.7 %, there is no deep reason why the 



finite parts in (5.1) and (^2) should be equal, as figures may vary by just changing the values of the different 
cutoffs. In a (fermionless) QCD context \x can arguably be related to the momentum scale parameter Aqcd (see 
e.g. Q and B). Results for that model follow by just bringing in an obvious (N% — l)-factor from the SU(iV c ) 
degrees of freedom; thus, for N c = 3 the cutoff-independent part of the energy becomes ~ (When fermions 
are assumed to be massless, neglecting them is not too bad an approximation, since results in the second work 
of refs.|2(J showed their contribution as being one order of magnitude smaller). 

Back to the e.m. case, the same confining set-up — without external modes — in cubic cavities (see |30| ||) 
yields a Casimir energy that amounts to 0.0916/L, where L is the edge length. The resemblance among this 



number and the finite parts of (5.1) , ( p.2| ) happens to be striking, although in view of the different details in the 



schemes leading to their derivation one ought to be cautious before taking this point any further. Joining (3.32) 



and Cm 1 ' 11 ' 7 * ( s ) i which is equal to ( |4.8| ), we shall find the net Casimir energy for a scalar field filling the whole 
space and satisfying Dirichlet b.c. on the spherical surface. Whithout having to apply (|T^), one gets a finite and 
scale- independent result, which reads 

Eg* 1, v = - 0.00282 (5.3) 

This finiteness is due to the cancellation of both poles at s = — 1 when adding up internal and external contri- 
butions, and may be put down to the self-erasing of curvature-dependent infinities with same size but opposite 
sign on each side of the surface. ( |5.3| ) coincides with the energy value which would yield the force (3.24) of ref. 
]To| | for the same physical situation, derived from a Green function approach. 

Finally we consider an e.m. field in the whole space with the sphere acting as a neutral and perfectly 



conducting boundary, which corresponds to the sum of the four results ( 3.31 ), ( 3.42 ), (4.8), ( 4.14 ). On addition, 
we realize that the poles cancel within the 2?-pair and within the 7?.-pair separately, rendering the PP prescription 



in (3.2) unnecessary as we are left with a finite and scale-independent value, namely 



i7ie.ni. 



— 0.09235 
2a 



(5.4) 



that coincides with the celebrated figure of Jij], |lj, |T^] . 

Let's take another glance at all the results for Eq prior to the application of PP. The outcome is summarized 
in table |l|. This way one easily sees that the PP prescription is redundant when any internal-external pair is 



D = 3 


Pirichlct 


TvLobin 


(a(3) = 1/2) 




1=0 

(only region I contributes) 


i r 7t i 

2a L 12 J 


1 IT 

2a 24 


{I > 1} region I 


1 

2a 


315 \ 
+ 


7TT + ln ( a M)) 
0.27069] 




1 

2a 


45 \ 


iTT + ln ( a M)) 
-0.10285] 




{I > 1} region II 


1 

la 


2 

_ 315 


0.00326] 


) 


1 

2a 


2 

_ 45 


(iTT + ln ( a ^) 
-0.07223] 


) 



Table 1: Zero-point energy decomposition in terms of scalar fields satisfying Dirichlet and Robin b.c. on a spherical 
surface of radius a in D — 3. 



added up, which may be envisaged as the above commented surface divergence cancellation. An analysis of 
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heat-kernel coefficients allows one to realize that these singularities are odd — and therefore of opposite sign on 
each side — when the space dimension is an odd number. In zeta regularization this is no longer so for curved 
surfaces in even D (see H [32| — notice also that the infinities in (To) confirm this observation). This fact has to 
be faced, in particular, in D = 2. 



5.2 D = 2 

Although we have argued that the e.m. problem in D = 2 reduces to a Neumann field, for completeness the 
analogous figures associated to a Dirichlet field were also calculated, and all the values obtained have been listed 
in table 0. Adding up all the Neumann parts, which are in the second column, one finds 



D = 2 


Pirichlet 


A/eumann (a(2) = 0) 


(regions I + II) 


£ [-0.02802] 


£ [-0.50704] 


{I > 1} region I 


i 

2a 


+0.02436] 


) 


i 

2a 


^ (^ + huV)) 

+0.17883] 




{I > 1} region 11 


J_ 

2a 


+0.00501] 




1 

2a 


-0.03804] 


) 



Table 2: Zero-point energy decomposition in scalar fields under Dirichlet and Neumann conditions on a circular line of 
radius a in D = 2 



a 



-JL]n( 0/i )_ 0.18312 



(5.5) 



This minus sign makes us think of the attractive force obtained for a cylinder in ]54[ | (although that refers to 
three-dimensional space, the symmetry is the same). Not even the sum of inner and outer parts cancels all the 
singularities. This is a consequence of using zeta regularization with a curved boundary in even space dimension 

ill- 

Doing the same for a scalar Dirichlet field we have 



£-al. = 1 

a 



— l-]n(a(i) + 0.00068 
128 



(5.6) 



The coefficient of the logarithmic part agrees with the first ref. of j5l| for the same set-up. Since the regularization 
method in that paper was frequency cutoff, this coefficient is all that one should expect to coincide. About possible 
ambiguities coming from the use of different regularization schemes, see e.g. the comments in [ |53"[ . 

A comparison with the results in ref. is now in order. The contributions from the I = mode in the 
Dirichlet (scalar) and Neumann (e.m.) cases agree with formulas (A13) and (3.5) — respectively — in that work. 
However, after performing the sum for infinite values of I, a divergent part — — formula ( |5.5| ) prior to 

prescription — comes into existence through the pole of Ch{z, 1) = Cr( z ) when z equals one (i.e. what we have 



called, below eq. (|3.24| ), a b-type singularity). Now, this specific piece has survived the infinity cancellations which 
take place when adding all the internal and external parts. Observing the parity of all our J7n(s)'s when going 
from internal to external parts, it is not difficult to realize that this happens for even D > 2. 

These poles were also detected in sect. Ill of ref. [ ^0| after using dimensional regularization and employing 
Riemann zeta functions in the last stage of the calculation. Within our zeta-regularization context, this divergence 
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would invalidate the alleged reliability of the finite estimates in ref. |5?J for the I ^ mode contribution (formulas 
(A16) and (3.9) of that paper, for scalar and e.m. cases, respectively). The existence of a singularity was there 
acknowledged, but it was argued that it could be a 'spurious' one. In our own framework, that speculation seems 
to be more questionable. By way of connecting results, we show in app. C how to reobtain those estimates by 
performing what might be called deliberately 'naive' zeta regularization. 

5.3 Ending comments 

Canonical field quantization leads to operators — e.g. the Hamiltonian — with ill-defined vacuum expectation 
values. Such troubles, due to quantum fluctuations, are often suppressed by decree removing them when no 
observable effects are expected. The picture is different in the presence of external sources or boundaries that, 
by breaking symmetries, render fluctuations observable. This explains the longevity of Casimir's general concept 
of vacuum energy, according to which the physical vacuum of quantum fields must be determined including their 
constraints. 

Full vacuum energies may contain nonobvious infinite pieces originated in boundary surface tensions, but our 
attention here has focused on the parts responsible for Casimir forces, i.e. those containing dependence on the 
relevant space parameter, which may be finite in spite of a global energy singularity. That is the way in which 
'finiteness' has to be qualified [Isfl . 

Viewed as an evaluation method for the Casimir effect, zeta function regularization had been relatively suc- 
cessful up to now but it was leaving some vague aftertaste of scepticism insofar as it was applied to comparatively 
simple problems: parallel plates, hypercubes, torii, hyperspheres as framework spaces (not as boundaries), i.e. 
situations where eigenvalues are at worst polynomials in the quantum numbers. Circular or spherical boundaries 
are beyond this realm. The unified approach of the present work has quickly enabled us to recover 

• three remarkable Casimir energy results involving a sphere: 

1. e.m. field inside: second row of table [|, which gives the coefficient for the logarithmic term in formula 
( |5.l| ) as in p0[ . The closeness of the remaining pieces may be in principle fortuitous since they come 
from different regularizations (but is numerically convenient for comparing scales). That is why we 
should regard our finite part as a new zeta-regularization result. 

2. scalar Dirichlet field inside and outside: first column in table [j], producing (5^3) like in ref. [jTo|| . 



3. e.m. field inside and outside: second plus third rows in table [j], that yield (5.4) coinciding with refs. 
and one concerning a circle: 



1. scalar Dirichlet field inside and outside: first column in table [| giving (5J3) whose logarithm coefficent 



is the same as in [^TJ . The rest is a new result of zeta function regularization. 

Furthermore, combining figures other findings emerge, e.g. the sum of all the contributions for the D = 3 
Robin scalar field — second column in table [l] — gives Eq 3,1 ' n = —0.02209—, of a larger order of magnitude than 
and opposite sign to its Dirichlet counterpart (|5.4|). Moreover, the zeta- function regularized version of the e.m. 



Casimir energy inside and outside a circle (5.5) is, to our knowledge, another unreleased result. 
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We hope that the new answers provided — together with the recovery of figures originally obtained after 
considerable effort — by this single zeta function strike can give the reader reasons to think that the versatility 
and scope of this technique may be somewhat wider. 



A Appendix: external modes 
A.l External Dirichlet modes 



Here we give a detailed proof of expression ( |4.2| ) . Our starting point is the defining series of the zeta function for 
a scalar field existing between two spheres. The inner one has radius a and is assumed to be physical. The outer 
one has radius R and it is introduced for the sake of convenience; eventually we shall take the R — > oo limit. 
Taking into account that v{D = 3,1) = I + |, initially we take solutions with radial part Ah^\Xr) + Bh^\Xr), 
where the hi's are spherical Hankel functions and A, B are coefficients whose relative value — with respect to each 
other — will be determined by imposing our Dirichlet conditions on both surfaces, i.e. at r = a and r = R. Doing 
so, we are left with an homogeneous linear system for (A, B) and, by requiring its compatibility, the equation 

f a>R (X) = /i ; (1) (Aa)/ij 2) (Ai?) - h^iXR^iXa) = 

follows. Therefore we study the 'partial- wave' zeta function £„(z) = — — , where A„ jP is the p-th zero of the 
function f a ,R(X). 

We tread here along the same lane that was open in , which runs through 

&(*) = 2^ / dX ^^MAflCA)). (A.l) 



In expression (A.l) the integration contour winds counterclockwise round the zeros of f a ,R, which are real. As it 



stands, this gives a representation for valid if Re z > 1. We slightly modify our representation to 

C^) = ^/<*AA— ^(l^^/^A)), (A.2) 

which ensures that the argument of the logarithm goes to 1 when A approaches 0; this will be seen to be useful 
in the sequel. 

Now we deform the contour in such a way that it is made up of three parts: a straight line from +ioo to ie, 
an arch of radius e connecting ie to —ie, and a straight line from — ie to — ioo. The courteous reader may care to 
see in [|6j that 

/«,ji(A) = e lXa - lXR S n {~iXa) S n (iXR) - e - lXa+lXR S n (-iXR) S n (iXa) , (A.3) 
fc=0 v ; 

With this in mind, let us try to find an analytic continuation to the contribution from the upper straight line. 
We immediately see that the dominant contribution in the argument of the logarithm comes from the e tXa - lXR 
term. So, we separate this factor in the logarithm and apply the property that the logarithm of a product equals 
the sum of the logarithms of its factors. This leads to 

,1-z 



^+W = -7r-7-rUi?-«) 6 



dp In 



pz+l 



1 va n R r ' 



2mZ+1 1 1 (A.4) 



pR 2v 



- (s n {pa)S n (- P R) - e- 2 ^ R -^S n (pR)S n (-pa)) 
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This expression explicitly defines an analytic function for Re z > 0. There is an equivalent representation for 
£„._ (the contribution from the lower straight line). It reads 



dp 

z+l 



hi 



P 



1 va n R 



zv 



2tt 



[R - a) ■ 



A-z 



1 



- (s n {pa)S n {- P R) - e- 2 P( R - a 1S n {pR)S n {-paj) 



(A.5) 



_p R? v - a 2 " 

The complete Q u is given by the addition of these contributions plus the integration along the small arch of radius 
e. In any case we have now an explicit continuation valid for Re z > 0. If we restrict to the domain < Re z < 1, 
and take the limit e — > we end up with 



C, v (z) = - sin (-2 



dp 



In 



1 i/n n F? n / \ 

- p R 2 V _ a 2u [S n {pa)S n {-pR) e-^ R -^S n ( P R)S n (-pa)) 



(A.6) 



By now it may seem unpleasant that the large- R limit is ill-defined. When we go on with the continuation process 
towards domains in the complex set with negative real part it will be apparent that in those regions one may 
perform this limit, which will allow us to extract physical results. 

In order to go on we have to add and subtract a function that apes the asymptotic behaviour of the integrand 
for large A and which may be analytically integrated, this is the general procedure which was set forth in |50| pH] . 
Nevertheless, before we do this, we proceed to further simplify our expressions. First of all, we separate the 
integral into three parts: 

dp 



I" JL P 



z+l 



hi 



z 



dp 

,2 + 1 



hi 



S n (pa)S n (-pR) - e- 2p{ - R -^S n {pR)S n {-pa) 
1 va n R n 



k Jl 9 
L dp 



,2;/ 



Z 



k Jq 9 



z + l 



In 



pR? v 



1 va n R r ' 



p R? v - a 2 " 



hi 



S n {pa)S n {-pR) - e- 2pi - R -^S n {pR)S n {~pa) 



(A.7) 



where L is any positive number. The second integration is trivial and gives a meromorphic function with a unique 
pole at z = 0. The third integration directly defines an analytic function in Re z < 1. It is the first integration 
which calls for the special treatment that has been overviewed above. So, we write 

GaiO) = ^sinQzJ 
-S n {pa)S n {-pR)-e- 2 P( R -^S n {pR)S n (-pa)A . .. , „A ( A " 



l 9 



dp 

,z + l 



In 



A n (pa,pR) 

where A n (pa, pR) is a function which shares the asymptotic behaviour for large p with 

S n (pa)S n (-pR) - e- 2 ^ R - a ^S n (pR)S n (-pa) 

in such a way that 



In (A n (pa,pR)) 



dp 

z+l 



In 



S n (pa)S n (-pR) - e- 2 P( R -^S n (pR)S n (-pa) 



(A.9) 



l P^ X L A n (pa,pR) 
defines an analytic function for — z Q < Re z (where z Q is a positive number which can be made as large as we 
wish by dint of complicating A n ), and 

dp 



J L -^t^Mp^pR) 



(A.10) 
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is easily computed, giving an explicit meromorphic function. Once this is achieved we have our much coveted 
analytic continuation of Q v for Re z > —z a . In this setting, we may restrict the function Q u to — z Q < Re z < 0; 
in this domain we have 



Uz) = C,i( z ) + J sin (f*) / ^^(s n (pa)S n (-pR)-e- 2 ^ R -^S n (pR)S n (-pa)) 
as happens to cancel one contribution from £„ 3 . To sum up, we write (— z a < z < 0) 



(A.ll) 



— sin — z 



— sin — z 



dp 
L dp 



In 



o P° 



S n (pa)S n (-pR) - e-M R - a )S n (pR)S n (-pa) 
A n {pa,pR) 

-In (s n {pa)S n {-pR) - e- 2 ^ R -^S n {pR)S n {-pa)} 

Z /7T \ f°° dp , A . 

-sm^-zj —^hxA n {pa,pR) 



+Analytic continuation of 
We have specified that we should calculate the analytic continuation in the last term, as the integral 



z fir \ f°° dp . 

7T V 2 Z ) J L ' ' 



(A.12) 



(A.13) 



only exists for Re z > — 1. Now, if A n (pa, pR) also satisfies that expression ( A.13 ) exists in the domain < 
Re z < 1 and its analytic continuation may be easily computed in the particular case that L — 0, then expression 



(A.12) may be simplified to 



In 



S n (pa)S n (-pR) - e-^ R - a )S n {pR)S n {-pa) 



+Analytic continuation of 



A„(pa,pi?) 
f Sin (f Z )/ ^^Mpa,pR) 



(A.14) 



The sagacious reader will immediately recognize that this expression admits a simple large- R limit, which may 
be written as 



■ sin 



dp 



In 



o 



Snjpa) 
An(pa) 



+Analytic continuation of 



z . fit 
— sin — ; 

7T V2 



P 



^ In A n (pa) 



(A.15) 



where A n (pa) is supposed to have the asymptotic behaviour of S n (pa) to some order. Note that this limit has 
been made possible when we have passed to regions with negative real part. To finish this digression, we cast 



expression (A.15) into a more standard form: 

Uz) = ~sm(lz)J o J^ln 



^e" a K n+h {pa) 
A n (pa) 



+Analytic continuation of 



~z ■ 






— sin 




f 


7T 




Jo 



P 



(A.16) 



If we only need a representation valid in — 1 < Re z < we need not include any A n at all and the result is 



«')-;-(5')f^ h (^*. 



Spa) 



(A.17) 



which, when rescaled to a = 1, and for n + 1/2 = ^(Z? = 3, n), gives (4.2). 
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A. 2 External Robin modes 



We explain now how to arrive at the expression for Robin boundary condition in the external region. 

In this derivation we will keep in mind what we did for the modes satisfying Dirichlct boundary conditions. 
So we shall also consider the existence of a large spherical shell of radius R which will be taken to oo once we 
have performed the analytic continuation to domains with Re z < 0. 

Note that we will be assuming here that the conditions that the modes satisfy on the outer sphere are also 
of the same Robin type as those which are imposed on the physical one. This is equivalent to saying that the 
space is confined by a large sphere which is a perfect conductor. In any case this particular choice of boundary 
conditions on the outer sphere has no physical relevance in the sense that when one has performed the analytic 
continuation and R — > oo, the results do not depend on it, for instance, we might as well decide to impose 



Dirichlet boundary conditions for both TE and TM modes on this unphysical sphere, and the final result (A. 25) 
would not be changed at all. 

Setting Robin conditions at r = a and r = R on the radial part of the spherical wave solution, one is led to 
consider the A„(3 n )P 's which are zeros of the function 



f a , R (X) = a^iXaWiXR) - a} lJ (XR)a^ J (Ao), 



r (2) 



(1), 



r (2). 



(A.18) 



where o\ l \z) = hn\z) + z4-hn\z) and u (3, 1) = I + i. In more detail one has 



r(2) 



i n+2 e~ lz 



-S n (iz) + S n (iz] 



As we did in the case of Dirichlet boundary conditions we first confine our efforts to the study of C, u (z) = 

E- — . which is again given by (A.l), but now with the f aR (X) in eq. (A.18). The integration contour winds 
\z 

p v >p 

counterclockwise round the poles of f a ,R, which are real (this property is easily drawn from the fact that these 
poles are eigenvalues of a self-adjoint operator). Our complex integral is a valid representation for if Re z > 1. 
We slightly modify it by inserting a harmless factor in the argument of the logarithm, which is again contrived 
so as to make the log an infinitesimal quantity when A approaches the origin of the complex plane 

uXa n+l R n+1 



z 

2~7ri 



dX x- 



ln 



(A.19) 



X n{n + l)i(R 2l/ - a 2 ")' 
this is valid if n > 0. 

Now we deform the contour the same way as we did before: a straight line from +ioo to ie, an arch of radius 
e connecting ie to — ie, and a straight line from — ie to — zoo. Now the procedure closely follows what we did in 
the case of TE modes. We find three contributions: , Cu,- an d the contribution from the arch. Both Cv,+i 
Cw- are given by 



U+(z) = ~ 

CvAz) = 



2tt 1 



2tt 1 - z 



-(R-a) 



-z-l 



(R-a) 



2tt 

:i z+1 
2tt 



dp 



hi 



up 



(I 



n+1 



dp 



In 



n(n + 1) R 2v - a 2 

,n+l j^n+l 



;e-rt R -^f a , R (ip) 



up 



n(n + 1) R 2v - a 2v 



e-^-o) f a , R {ip) 



(A.20) 
(A.21) 



It is immediate that we have realized an analytic continuation for Re z > 0. 



22 



If we restrict the domain to < Re z < 1, and take the limit e^Owe end up with (note that in this limit 
the contribution from the arch is vanishingly small) 



£,(*) = J sin (| 
For the sake of clarity we give here 



o P 



dp 

z + l 



111 



n+1 pn+1 



-vpa n+L R 



n(n + 1){R 2 » - a 2 ") 



fa, R Hp)e- piR - a) 



(A.22) 



- {R - a) faA-w) = (s n (pa)-S' n (pa))(s n (- P R)-S' n (-pR)) 

(j5 n (pR) - S n {pR)^ (S n (-pa) - S n (-pa) 



-2(R-a) 



(A.23) 



Now the reader should have no difficulty to apply the same procedure that we have explained in the case of 
external TE modes. In particular, we will have that for any z there is a proper function A n (pa, pR) such that 
the analytic continuation of Q v for — z Q < Re z < is given by 



Analytic continuation of 
dp 



+ 



Z . /7T \ 

— sin — z 
7T \2 J 



o P 



z+l 



In 



^ Sin (f Z )/ -^^A n {pa lP R) 
-e-^f a ,R(-ip) 



A n (pa, pR) 



(A.24) 



The choice of function A n is not uniquely determined. By simple inspection one sees that it would be enough 
a function such that for large p, A n (pa, pR) has the same (algebraic) asymptotic behaviour (to some convenient 
order) as the product 

(S n (pa) - S' n (pa)j (s n (-pR) - S' n (-pR)j . 

A proper choice might be of the form A n {pa, pR) — An (pa) An\pR) , where An\pa) is an algebraic function 
that asymptotically behaves like 



S n (pa) - S n (pa) 



(2) 

and An (pR) does the same job for 



S n (-pR) - S n (-pR). 

This would suffice for our purposes. As liniR^oo A n ' (pR) = 1, which is easily concluded from its asymptotic 
behaviour, we see that in the domain which we are considering it is possible the R — > oo limit. This limit leads 
to 



Cv(z) = Analytic continuation of 

dp 



Z . /"IT 

— sin — , 

7T V2 



dp 



■ sin — z 



o P 



In 



'o P 
S n {pa) - S' n (pa) 
A n 1] {pa) 



z + 



T In A£\pa) 



If only a continuation valid in — 1 < Re z < is needed, then we simply have 



S n (pa) - S n (pa) 



and, for a = 1, 



Uz) = ^singz) 



o P 



dp 

z+l 



111 



e z -zK v (z) - -K v {z) 



(A.25) 



(A.26) 



(A.27) 



which gives (4.11) after noting that a(D = 3) = 1/2. 
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B Appendix: heat-kernel series approximation to the complete spec- 
tral zeta function 



For the Laplacian □ (or any such operator satisfying some suitable requirements) the spectral zeta function 
(n{z) = Tr U~ z is related to the heat kernel Yb(t) = Tr e~ tD through the well-known Mellin transform 

1 



r ( z ) Jo 

In turn, the heat kernel has small-f asymptotic expansion 

1 



(A-Kt)°/ 2 



dtt z - L Y n (t). 



k>0 



(B.l) 



(B.2) 



We split the integration domain into the intervals [Q, r] and [r, oo), where r is such that (B.2) holds for t < r. 
Thus, one can use this series in the first part and integrate, finding 



1 



T(z) 



1 



E 



z+k/2-D/2 



(4nt) D / 2 ^ Q z + k /2-D/2 



/CO 
dtt z - l Y n {t) 



(B.3) 



In our case, where the boundary is a sphere, the values of the &fc/2's can be read off from |L8| (up to k — 20 if 
necessary!). Since □ gives squares of e.m. modes, symbolically □ = M 2 and (m( s ) = (n{s/2), which shall be 
studied around s = — 1. Thus, separating the singular part at this point we can write 



Cm(s) 



s+1 



+ p{t)+c(t) + 0(s + 1), 



(B.4) 



where 



2b, 



(D+l)/2 



p(r) 



c(t) 



(47r)^r(_i; 
l 

i 



(fe-Z3-l)/2 



W)/ 2 (lnr-V^(-l/2)) + 2 ^ h _ D _\ 



fc>0 



(B.5) 



From heat-kernel properties, one also knows that for large t 



Y n (t) 



-Ant 



(B.6) 



where Ao stands for the smallest eigenvalue of □. We wish to use this behaviour in order to approximate the 
value of c(t). 

As an example, we take D — 3 and a scalar Dirichlet field inside the sphere. From the bn listed in e.g. |lS| , 
we find 

as we already knew (( 3.3 2| ) ) . Ao will be the square of the smallest nonvanishing zero of Jx/2{x) which is x = ir; 
therefore Ao = ir 2 . Now, a delicate question arises about how to choose r large enough so that the replacement 
(B.6) be sensible while maintaining the validity of the small-i expansion (B^) for t € [0,t]. The answer will 
necessarily be a compromise between both requirements. Looking at the integral c(r), which we hope to keep 
small, one considers the region near t = r and realizes that it would be desirable to have Aqt 3> 1. Thus 
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largest k 


approx. p(l/2) 


5 


0.009230 


10 


0.009331 


15 


0.009333 



Table 3: Approximations to p(l/2), depending on the largest k included in the sum. 



one arrives at — -C t < 1. Since Ao = tt 2 , we deem r = 1/2 as a fairly adequate choice. Then we seek an 
approximation to c(l/2), which is shown by means of table BL On the other hand, the approximate value of the 



integral is, by (B.6), c(l/2) ~ —0.0004. As a result, the finite part of C>f( s ) ls 



p(l/2) + c(l/2) ~ 0.0093 - 0.0004 = 0.0089 



(B. 



in good agreement with the finite part of (3.32) 



C Appendix: 'naive' zeta-function regularization 

Let's calculate again the e.m. case in D = 2, but making now the addition of internal and external modes before 



the angular momentum summation. After adding up ( 3.35| ) to (4-11) for a = (i.e. purely Neumann b.c.) and 
rescaling x — > vx, we find 

[£*(.) + C?'*(-)] a= o 

dxx'^ 1 \Yi[-2vxI' v (yx)K' v {vx)\ . 



S . TVS 

— sin — v 

7T 2 



(C.l) 



Making use of the u.a.e.'s for V v (yx) and K' v {yx), taking advantage of previous steps and of parity reasoning 

1 



In [-2vxl' v (yx)K' v {vx)\ - In 



xt(x) 



E 



2Z4f(*(x)) 



a=0 



(C.2) 



where t(x) is the one in (3.14) and the Wj^'s are given in (3.37). Subtracting up to n = N terms from the 
integrand, one may write 



„\r , > s . its 

IT A 



dx x 



-i„" + ia » 1/a + v 



(C.3) 



By ( |3 .11 ), we know that the above integral amounts to 2 • - — : — — The J^^s are already specified by the 

sin 2 

integrals ( 3.4C ). Sff(s;iy) denotes the corresponding residual integral, in the style of ( |3.3S ), containing the 
difference between the exact integrand and its iV-term approximation. 

Taking into account the D = 2 degeneracies and noting that f(2,l) = I, the complete zeta function is put as 



C#,(*)=« s 



1=1 



The I = contribution to the Casimir energy is therefore 

itfc-D -£[<*«(.) +<?■«<.) 



ct=0, s=-l 



-0.25352 

a 



(C.4) 



(C.5) 
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(formula (4.16)). This agrees with (3.5) in ref. |37j]. On the other hand, the set of I ^ modes yields a 
contribution which amounts to the s — > — 1 limit of the following quantity: 

N 



;=i 



[ 2V oo 

-Cb(«) + - sin ^ E 2 C«( 2 " + s ) # («) + E S n(^ ' 



(C.6) 



At s 



-1 the Z-sum causes no problem, as it is finite and relatively small. However, the n-sum contains a 



divergence in its n = 1-term, coming from the pole of when its argument equals one. Bearing this in mind, 



we find the residue of ( |C.6| ) at s = —1 to be —J 2 ' {— 1) = , as given by (5.5) (up this point, we have 

just been doing the same calculation but re-grouping things in a convenient way). Now, if we deliberately ignore 
this singularity and throw away as 'corrections' the I— and n— sums in (|C.6|), we simply get 



2a ZAa 



(C.7) 



which is precisely the 'LT' ('leading term') contribution of the I ^ modes to the Casimir effect given in eq. 
(3.9) of ref. J37J, i.e. it is the part which comes from keeping just the first term on the r.h.s. of (C.2). Of course, 



viewed from inside our zeta-regularization context, such an approximation in unjustifiable, since the whole value 
has a divergence attached which we dare not call 'spurious'. 



Doing the same with the scalar Dirichlet modes, one easily sees that the I = part is one half of (4.10), in 
agreement with result (A13) of ref. J37J, and the I ^ piece analogous to the above mentioned approximation is 



Za 24a 



(C.8) 



which amounts to the 'LT' result (A16) in ref. p7[ and whose validity rests on the same flimsy assumptions. 
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